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A multimode optical fiber with a truncated transverse cross section acts as a powerful versatile 
support to investigate the wave features of complex ray dynamics. In this paper, we concentrate on 
the case of a geometry inducing a mixed dynamics. We highlight the properties of regular modes 
associated to stable periodic orbits such as an enhanced intensity localization and report unexpected 
attributes such as the statistics of the Inverse Participation Ratio that present features analogous 
to those of Anderson localized modes. Our study is supported by both numerical and experimental 
results. 
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I. INTRODUCTION 

Chaotic systems are more and more used in applica¬ 
tions. These systems, whose geometry induces complex¬ 
ity through a chaotic ray dynamics, bring a new kind 
of solutions for practical purposes. Recent developments 
in the field of electromagnetic compatibility show that 
chaotic geometries promote a very efficient mode mix¬ 
ing in electromagnetic chambers jl, 2]. In optics, dielec¬ 
tric microcavities stand as the typical example of the ex¬ 
ploitation of chaos as an improvement of existing devices 
[3|. Indeed, chaotic microcavity lasers have shown to al¬ 
low low threshold and highly directional outputs i@ 
and are a very useful solution for integrated optics. Re¬ 
cently, microcavity-based gyroscopes in which a chaotic 
cavity enhances the rotation sensitivity ji[ have been in¬ 
vestigated. Before that, a chaotic double-clad fiber am¬ 
plifier had been proven to be an optimized solution for a 
uniform pump absorption along the total amplifier length 

0 . 

Instead of considering the fully chaotic regime, one can 
consider a so-called mixed system. In this case, the co¬ 
existence of regular and chaotic motions leads to an en¬ 
larged wealth of the dynamics. This duality gives for 
instance rise to a wave feature known as dynamical tun¬ 
neling, which allows a passage between the chaotic and 
regular regions by coupling one to another [8]. This effect 
has been used in segmented waveguides where tunneling 
between stability islands ensures a non-diffractive beam 
regime [9], as well as in deformed microdisks in which 
light is coupled to a bus waveguide through the tunnel¬ 
ing effect acting as a resonant dynamical filter 0- 

A highly multimode optical fiber whose transverse 
cross-section presents a truncated-circle shape allows to 
explore different regimes of the dynamics by changing the 
length of the truncated diameter (e.g. changing d in Fig. 
DJb)). Over the past few years we studied experimen¬ 
tally manifestations of ergodicity as well as deviations 
from this universal behavior due to scar modes in pas¬ 
sive and active chaotic optical fibers HM3- We keep 
on exploiting the potential of multimode fibers, but now 
concentrating studies on wave features of mixed systems. 


The latter present both regular and chaotic dynamics 
in the ray limit, as exemplified by the dynamics of the 
so-called Mushroom billiard [Q~5]. In this paper, we fo¬ 
cus on the properties of the modes of a mixed dynamics 
based multimode optical fiber, and in particular modes 
that are associated to the regular part of the dynamics. 
Unlike ergodic modes whose statistics follow a universal 
behavior, the regular modes exhibit some strong devia¬ 
tions that are unexpectedly analogous to those induced 
by localized modes in disordered systems 0 - 

In section im we present our experimental system. We 
also present two usual manners of distinguishing regu¬ 
lar, chaotic and mixed systems, on the one hand through 
their Poincare Surface of Section and on the other hand 
through their spectral statistics. In section m we de¬ 
velop the analysis of the spatial properties of the regular 
modes of the fiber using the analogy with the gaussian 
modes of an optical resonator. In section [IV] we pro¬ 
pose an original point of view by interpreting our results 
through tools usually used in disordered systems exhibit¬ 
ing (Anderson) localized modes. In section |Vj we present 
an analysis of the properties of the regular modes of the 
optical fiber based on both experimental and numerical 
results. Then, in section ED we conclude and propose, 
in an open discussion, some potential applications. 


II. THE OPTICAL FIBER AS A VERSATILE 
OPTICAL ANALOGUE OF A DYNAMICAL 
SYSTEM 

Our experimental system is a non-standard silica mul¬ 
timode optical fiber, whose cross section is a truncated 
circle. The diameter of the core of the fiber is 2 a = 
125/im, and the truncated diameter d = 7 x a, with 
7 G ]0; 2]. The core is surrounded by a silicon cladding 
of diameter 250/im [Fig. [Tj . The optical indexes are 
n co = 1.458 and n c \ = 1.41 for the core and the cladding 
respectively for the vacuum wavelength Ao = 632nm. As 
the wavelength is small compared with the characteris¬ 
tic size a of our system, the semi-classical approxima¬ 
tion is valid, and one can consider both the classical de- 
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FIG. 1 . Presentation of the experimental system, (a) picture 
of the cross section of the optical fiber for 7 = 0.95 taken with 
a standard microscope in the transmission mode, (b) scheme 
of the fiber with 7 < 1. 

scription of ray dynamics as well as a modal description 
of propagative waves. In this regime, the cross-section 
of the optical fiber is equivalent to a two-dimensional 
(2D) bounded cavity which is a well-known paradigm of 
a hamiltonian dynamical system (2D billiard) to study 
wave chaos 0- Indeed, the longitudinal evolution of a 
ray along the fiber is equivalent to the time evolution of 
a trajectory in a billiard. Thus, in the following, we will 
refer to the transverse cross section of the optical fiber 
as a “cavity”. 

The versatility of our system relies on its ability to ex¬ 
plore different types of dynamics from regular to chaotic 
by varying the factor 7. The qualification of the dynam¬ 
ics is commonly studied through a stroboscopic phase- 
space representation, well-known as the Poincare Surface 
of Section (SOS). The SOS is a reduction of the ray dy¬ 
namics to a 2-dimensional mapping. It consists in plot¬ 
ting at each impact of the ray on the boundary the curvi¬ 
linear abscissa 8 measured through the angle 0 = s/a and 
the sine of the incidence angle, a [Fig. Efa)]. As shown 
in Fig. Hc-e) and as pointed out by Ree [18|, the SOS 
presents the characteristic behaviors related to different 
types of dynamics while varying 7. The nature of the 
classical dynamics in a D-shaped billiard as been proved 
to depend on the stability of the shortest periodic or¬ 
bit (2-bounce periodic orbit (2-PO)) along the diameter 
d. When 0 < 7 < 1, the truncated diameter d is lower 
than the radius of the billiard. In such a configuration, 
the dynamics is mixed and the SOS presents a complex 
combination of stochastic and regular regions. In figure 
[Fig. Efc)] plotted for 7 = 0.95, regular islands coexist 
with diffuse points associated to the so-called chaotic sea 
0- Each regular island corresponds to quasi-periodic 
trajectory in the vicinity of a stable periodic trajectory 
at the center of the island. These stability islands mark 
the regular dynamics while the chaotic sea is a signature 
of chaotic trajectories. A large fraction of SOS is occu¬ 
pied by the main central resonance that corresponds to 
the 2-PO which is stable in this case. As d reaches the 
value of the radius a (7 = 1 , half-circle), the SOS tes¬ 
tifies of the regular nature of the dynamics [Fig. Efd)]. 


In this case, the angle of reflection a remains indefinitely 
the same for a given initial condition. Note that the 
SOS of a circular billiard (7 = 2 ) would be the same due 
to angle conservation. When d becomes greater than a 
(1 < 7 < 2 ), no more regular islands exist, and all the 
periodic orbits (PO) become unstable. Then, the SOS 
is densely covered by diffuse points and the dynamics is 
chaotic [Fig. [2](e)] . It is worth mentioning that the latter 
geometry has also been fully studied in optical fibers over 
the past few years 00 . In particular, we experimen¬ 
tally demonstrated that some specific scarred modes of a 
chaotic optical fiber - spatially localized along the least 
unstable periodic orbits - can be selectively enhanced 
through an optical amplification process 0 20 ]. More¬ 
over, as the pump absorption is clearly improved by the 
use of a chaotic double-clad fiber amplifier, the device has 
been proposed as an optimization of existing amplifiers 

UM- 

We now focus our study on the system with a mixed 
dynamics, and thus consider the case 7 = 0.95 which cor¬ 
responds to our fabricated fiber. The length of the fiber 
is 10 cm, which is greater than the Heisenberg length 
0 ensuring the validity of a modal description of light 
propagation. The number of modes of the fiber is given 
by the standard formula [ 22 | and is evaluated to approx- 
imatively N = 6000 at Ao- As shown on four exemples 
in the Figure 03 the modes exhibit some signatures of the 
underlying classical mixed dynamics (23[. The spatial 
distribution of intensity (near-field, NF) of the mode in 
Fig. [3](a) is located on the 2 -PO represented in pink in 
Fig. [ 2 ](b) . As mentioned before, this 2-PO corresponds 
to the central point of the major stability island in the 
SOS [Fig. EJc)]. Figure [3](b) shows the far-field (FF) of 
the same mode. The FF is the square modulus of the 
spatial Fourier transform of the field distribution. It in¬ 
dicates the direction R and modulus k of the transverse 
wave vectors. For Fig. [3](b) , two maxima of intensity are 
localized on the k x -direction of the 2 -PO at a distance 
k of the center of the Fourier space. This directivity 
in the FF evidences the signature of the underlying 2- 
PO (see Fig. [ 2 |). Strong correspondances between some 
particular modes and their associated classical trajecto¬ 
ries in the semiclassical regime are established through 
the Husimi representation, which is commonly used as 
a wave equivalent to the classical SOS 0 . In the case 
of 2D cavities, the Husimi function measures the nor¬ 
mal derivative of the eigenfunction on the boundaries of 
the cavity fl8l, I 25 I [27^ . Here, the Husimi function is cal¬ 
culated for the modes of the optical fiber and is repre¬ 
sented in the 2 D space (0,since). Here, 0 is the curvi¬ 
linear abscissa normalized to the radius as described in 
Figure [2I[a) ; sin a is the projection of the transverse wave 
vector R along the tangential direction to the boundary, 
normalized by the modulus k, = \R\ of the mode. Fig¬ 
ure He) presents the Husimi representation of the mode 
shown in Fig. [3](a) with the SOS of the 2 -PO super¬ 
imposed. The localized pattern of the intensity in the 
Husimi representation means that the mode is strongly 
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FIG. 2 . (a) Curvilinear abscissa s = aO and angle of reflection a used for the calculation of the SOS and Husimi representations, 
(b) Three examples of stable periodic orbits of the mixed billiard (7 = 0.95) : in pink the 2-PO, in blue the 5-PO and in green 
the 6 -PO.(c)-(e) Poincare Surface of Section representation for three different values of 7. (c) 7 = 0.95, mixed dynamics with 
a coexistence of regular and chaotic trajectories ; the pink star points out the 2 -bounce stable periodic orbit, the blue squares 
correspond to the central point of the stability islands of the 6 -bounce stable periodic orbit and the green triangles to the 
5-bounce stable periodic orbit, (d) 7 = 1 , regular dynamics for which one sees the conservation of the angle of reflection, 
whatever the initial condition is, (e) 7 = 2/3, chaotic dynamics, no stability regions, all the periodic orbits are unstable. 


confined along the x-direction. This confinement as well 
as the very good agreement between both this represen¬ 
tation and the SOS confirms that this mode builds on 
constructive interferences along the 2-PO. In the same 
way, Figures [3(d,e,f) and (g,h,i) show the NF, FF and 
Husimi representations of modes built in the vicinity of 
the 5-bounce [Fig. [2(b), blue line] and 6-bounce [Fig. 
mb), green line] stable periodic orbits respectively. As 
before, the agreement between the Husimi and SOS rep¬ 
resentation reflects the fact that the ray dynamics under¬ 
lies the modal behavior. Figures [3(j,k,l) present the NF, 
FF and Husimi representations of a chaotic mode of the 
fiber. The main difference is highlighted by the Husimi 
representation which is no longer confined in small areas 
delimited by the corresponding stability islands of the 
SOS, but spread over the chaotic regions. 

As conjectured by Percival [HI, the duality between 
chaotic and regular dynamics is also encountered in the 
spectral features. A simple way to quantify the degree 
of regularity of a 2D-cavity is given by the study of the 
statistics of energy level spacing 


sa 




n+1 


- E n 


N Z^n=1 - Ljn 


(i) 


where E n is the energy of mode n and N the total num¬ 
ber of modes. As described in [HI, HH, the probability 
distribution P (s a) in a mixed system is given by 


P(s A ) 


ds 2 A 


exp (— Wsa) erfc 




(1 — W) S A 


( 2 ) 

where 0 < W < 1 is the density of regular states. A for¬ 
mal evaluation of W is deduced directly from the SOS, 


by measuring the area of the regular region, weighted by 
the lengths of the associated trajectories. Figure 0] shows 
the statistics of a chaotic cavity (7 = 3/2) for symmetric 
[Fig. 0(a)] and antisymmetric [Fig. 0(b)] modes and a 
mixed cavity (7 = 0.95) for symmetric [Fig. 0(c)] and 
antisymmetric [Fig. 0(d)] modes. The gray scale his¬ 
tograms represent the energy level spacing distribution 
calculated for the modes of both cavities and the red line 
is a fit of the histograms by the function P(sa) giving 
an estimation of the parameter W. For the chaotic case, 
W = 0 indicates that the relative size of the regular re¬ 
gion is reduced to zero. This result is consistent with 
the typical SOS of a chaotic system [Fig. [2(e)]. Indeed, 
the inherent nature of a chaotic dynamics is to avoid any 
stable structure. On the contrary, for the mixed system, 
W ~ 0.5 corresponds to a SOS half-filled with the regular 
region. Note that this is qualitatively coherent with the 
SOS presented in Fig. He). The case W — 1 would cor¬ 
respond to regular systems 7 = 1 (half-circle) and 7 = 2 
(circle). 


III. REGULAR MODES AS MODES OF A 
PLANO-CONCAVE RESONATOR 

We now go one step further in the analogy between 
longitudinal optical ray propagation in the fiber and time 
evolution of the ray in the billiard (here the transverse 
cross section of the fiber). As evidenced in details by H. 
E. Turecci in [u| in dielectric microcavities, the modes 
localized on stable POs can be described in terms of gaus- 
sian optical theory. Using the standard physics of opti¬ 
cal resonators, we develop that correspondence for the 
2-PO regular modes as the analogous of the fundamental 
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FIG. 3. Near-field, Far-field and Husimi representation of 
(a,b,c) a regular mode of the 2-PO. Super-imposed in red is 
the SOS associated to the 2 -PO ; (d,e,f) a mode built on the 
5-PO. In red, the SOS of the 5-PO ; (g,h,i) a mode built on 
the 5-PO. In red the SOS of the 5-PO ; (j,k,l) a chaotic mode. 



FIG. 4. Probability distribution of the energy level spacing for 
the symmetric and anti-symmetric modes of (a,b) a chaotic 
cavity (7 = 3/2) and (c,d) a mixed cavity (7 = 0.95). 


longitudinal modes of a plano-concave optical resonator 
with dielectric boundaries where x acts as the longitudi¬ 
nal axis of propagation and y as the transverse direction. 
Thus, the optical resonator formalism allows us to ana¬ 
lyze both the stability of the cavity and the spatial field 
distribution of the modes of the fiber. 

First, the ray transfer matrix formalism enables to 
study the stability of the cavity by investigating the 
paraxial ray evolution in the vicinity of the 2-PO. Ray 
transfer matrix theory gives the state of output quanti¬ 
ties s out and a out in terms of input quantities Si n and c^ n 
using the optical characteristics of the resonator and its 
boundary conditions. Here, s refers, as described earlier, 
to the position of the beam on the resonator boundary 
and a to the angle of the beam compared with the normal 
on the boundary. In the paraxial approximation, output 
quantities are linearly dependent on input quantities : 


I Sout 

\(%OUt 



( 3 ) 


where M t is the ray transfer matrix. The stability of the 
resonator is given by the value of the trace of M t asso¬ 
ciated to a periodic sequence of a given ray. When rays 
are periodically refocused, the sequence is stable and ge¬ 
ometric characteristics of the resonator obey the relation: 


0 < 




( 4 ) 


where d = 7 x a is the resonator length and i?i, R 2 are 
the radii of curvature of the resonator boundaries. This 
stability condition is analogous to the stability analysis 
of dynamical systems based on the evolution of a small 
deviation in the vicinity of a PO [30j . The transfer matrix 
is then formally equivalent to the so-called monodromy 
matrix. 

In our case, R\ = a and R 2 = 00 because of the flat 
reflective boundary. From the relation one sees that 
the stability of the resonator depends on the value of d. 
Cases 7 = 1 and 2 are respectively associated to the hemi¬ 
spheric and confocal stable resonators. For 1 < 7 < 2, 
the resonator is unstable. For 7 < 1, the condition (j4]) 
is fulfilled so that the resonator is stable. The mixed 
fiber with 7 = 0.95 can thus be considered as a stable 
optical resonator whereas the chaotic fiber (7 = 3/2) 
corresponds to an unstable resonator [18|, [28|, [3l| . As a 
matter of fact, the nature of the dynamics of the under¬ 
lying billiard with a truncated circle shape is conditioned 
by the stability of the 2 -PO as mentioned before. 

In the following, we consider the case of the stable res¬ 
onator with 7 = 0.95, corresponding to the cross section 
of our optical fiber. We use the gaussian wave formalism 
in the (x, y )-plane in order to study the spatial distribu¬ 
tion of the modes built on the 2-PO. The spatial expan¬ 
sion of a beam traveling in the x-direction in the slowly 
varying envelope approximation reads: 


ip(x, y) = u(x, y ) x exp {-jux) 


( 5 ) 
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where u(x,y) is a slowly varying function describing the 
deformation of a plane wave due to the curvatures of the 
resonator walls and i/;{x,y) is a solution of the transverse 
Helmholtz scalar equation A'tp + k 2 ^ = 0. By using ex¬ 
pression (j5j) for ^){x,y), one gets the following Fresnel 
(or parabolic) form for the wave equation in the paraxial 
approximation: 


d 2 u d 2 u . du 

d^ + d^~ JK fa = 


(6) 


The solution u{x,y) of (| 6 |) — under the assumption of 
u{x, y) varying slowly enough with x so that d 2 u/d 2 x <C 
\2ndu/dx\ and under the condition of a matching be¬ 
tween the wave front and the curved and plane walls of 
the resonator - is given by [32j : 

u(x,y) = ^p-H m (V2-^—\ (7) 

w(x) \ w{x) J 

y exp - <*,„,) - P 2 } 

where w(x) is the beam radius, wo is the minimum beam 
diameter - the so-called waist diameter - for which the 
phase front is plane, R{x) the radius of curvature of 
the wavefront and H m is the Hermite function of order 
m. The solution u{x,y ) is consistant with the theory 
of gaussian beams along the ^/-direction [32|. The gaus- 
sian beam profile is characterized by the beam radius 
w(x). An initial gaussian beam of width wo at x = 0 
experiences a transverse expansion given by w 2 (x) = 

Wq ^1 + (x/x r ) 2 ^J during its propagation along the x-axis. 

x Y is the Rayleigh length, measuring the spatial coherence 
of the beam along the axis of propagation x. In our spe¬ 
cific case, the relevant wavelength is A^ = 2tt/k. Then, 
in our case, the Rayleigh length x Y and the beam waist 
wo are linked by the relation x Y = kwq/2. The beam 
propagating along the 2-PO also undergoes a phase shift 
which is twofold. First, the reflection on the dielectric 
core/cladding interface - governed by Fresnel reflections 
laws - implies a phase shift (j>* which depends on k. Sec¬ 
ond, the Gouy phase <p g [32, [33[ appears when the beam 
focuses. It results a complex total phase shift <fi t that 
reads : 



FIG. 5. Modes of the mixed fiber associated to different values 
of {m,p). (a)-(d) p = 4 and m — 0,1, 4, 7, (e)-(h) p — 12 and 
ra = 0,1,4, 7. 


Each mode of the resonator is then defined by the value 
of its transverse wave vector which depends on two 
integers, m and p running from 0 to maximum values 
satisfying < &o\/(nj?o - n^) 0 - In the optical res¬ 
onator analogy, p defines the number of nodes along the 
axial direction x , that is the order of the longitudinal 
mode along the 2-PO whereas m is the transverse mode 
number associated to the so-called high-order modes. 
The fundamental gaussian mode, that is the fundamen¬ 
tal mode for the transverse oscillations, corresponds to 
m = 0 and p — 0 . Some examples of the modes of the 
fiber/resonator with their corresponding value of (m,p) 
are given in Figure [5] 


IV. REGULAR MODES AS LOCALIZED 
MODES INDUCED BY THE MIXED DYNAMICS 


4>t=4>r + (f> g (8) 

= 2 arctan \J (n ^ 0 — n^) — 1 

+ (m + 1 ) arctan 

where ko = 27r/Ao is the modulus of the vacuum wavevec- 
tor. Stationary modes occur when the accumulated phase 
shift along a round-trip in the resonator is a multiple of 
27r. From 0 and (JEj), this condition of constructive in¬ 
terferences leads to: 



Km# x 2d = 2irp + 2<fir + 2{m + 1 ) arctan 
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(9) 


The ergodic modes of a chaotic system stand for the 
generic behavior. They are characterized by a gaussian 
statistic of the spatial field distribution. Nevertheless, 
some modes, namely seen as an “extra density (that) 
surrounds the region of the periodic orbit ” [34] are con¬ 
structed along unstable periodic orbits. They are usually 
named “scar modes”. Regular modes of a mixed system 
also result from constructive interferences although they 
take place along a stable periodic orbit. This intrinsic dif¬ 
ference in the nature of the underlying periodic orbit im¬ 
plies, among other features, a strong spatial localization 
of the field [Fig. [Sfa) and (e)]. This localization implies a 
substantial deviation in the statistics compared to the ho¬ 
mogeneous field distribution. The Inverse Participation 
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IPR 


FIG. 6 . Probability distribution of the Inverse Participation 
Ratio, (a) Chaotic cavity (7 = 3/2). The probability distri¬ 
bution is centered around the value IPR=3. (b) Mixed cavity 
(7 = 0.95) presenting some larger values of the IPR. The 
continuous line is a fit based on a nonlinear sigma model [35[. 


Ratio (IPR) and its statistics, which are commonly-used 
tools in the characterization of the spatial signatures of 
chaotic or disordered systems, notably allow to highlight 
this deviation. The IPR, that is the second order moment 
I 2 of the intensity, can also be defined as : 


(JJ s I(x,y)dS) 2 


where I(x,y) = \u(x,y)\ 2 is the field intensity, S the 
surface of the transverse cross section of the fiber and 
4 ff s I(x,y)dS = 1. The distribution of the IPR [Fig. 
[6] for the modes of the chaotic and mixed fibers ex¬ 
hibit some specific features. P{1 2 ) for the chaotic fiber 
(7 = 3/2) presents a peaked distribution around the 
value I 2 = 3 [HI, as expected by the universal prop¬ 
erties of ergodicity of the generic modes of a chaotic sys¬ 
tem [ljJ. Even if scar modes present some enhancement 
of intensity, this is marginal enough not to appear in the 
distribution of the IPR [Fig. [6(a)]. The Random Matrix 
Theory (RMT) predicts that there should be no fluctua¬ 
tions, but as the system is bounded, fluctuations appear 
[HI- On the contrary, the distribution of IPR calculated 
for the modes of the mixed fiber shows an asymmetric 
profile with enlarged high values up to 2 times the actual 
mean value (I 2 } — 4.5. This points out that a large num¬ 
ber of modes contributes to this deviation by present¬ 
ing highly localized intensities. In disordered systems, 
this behavior is currently associated to the presence of 
so-called “localized modes” some theoretical stud¬ 

ies based on the supersymmetry method |38j predict the 


asymmetry of the IPR distribution in the (Anderson) lo¬ 
calization regime. The distribution of the IPR is then 
expected to follow : 

P(I 2 ) = Cx yT X exp(-^/ 2 ) (11) 

for I 2 > (h) where ( I 2 ) is the mean-value of the IPR, 
C is a normalization constant and g the conductivity of 
the system that depends on the system size, ( I 2 ), and 
the mean-free-path l (HI- The high values of I 2 in Fig. 
[ 6 ]( b) present a very good agreement with eq. m with 
g = 2.02. This value is associated to a mean-free-path 
l = 16.8/im. The mean-free-path here is smaller than 
the typical size of the system whereas it is expected to 
approach the size of the system in the ergodic case. 

To investigate the origin of this deviation in the tail of 
the IPR distribution, we report in Figure [7] the value of 
the IPR of each individual mode as a function of^xa for 
the chaotic [Fig. [71(a)] and mixed [Fig. [71(b)] fibers. As 
expected for the chaotic case, the IPR tends to an asymp¬ 
totic value of 3 for most of the modes except for a small 
amount of them. The modes presenting larger values of 
the IPR [Fig. [7(a)] are some specific spatially localized 
modes built upon quasi-stable trajectories which follow 
neutral orbits. They are called bouncing ball modes, and 
take place along the full diameter of the fiber [see inset 
A in Fig. [7(a)]. They thus can be seen as a reminis¬ 
cence from the regular circular billiard and do not follow 
the standard gaussian statistics. Figure [7(b) presents 
the IPR as a function of k x a for the mixed fiber. In 
a surprising way, one observes a striking structuration 
of the values of the IPR for individual modes. More¬ 
over, one distinguishes some families of modes following 
a monotonous behavior. The upper red circles all cor¬ 
respond to the family of regular modes along the 2-PO. 
The high values of the IPR confirm a strong spatial lo¬ 
calization of intensity of these modes, also labelled as the 
fundamental gaussian modes uq^ p . The dots in shades of 
colors from red to green correspond to the IPR calcu¬ 
lated for the higher order modes u m ^ p with m varying 
from m m 1 to m = 9. The regularity of the IPR as 
a function of k x a is also observed for high er values o f 
m, the maximum being given for = &o\/(^co ~ n l\)> 
but they are not pointed out in the figure for the sake 
of clarity. The purple crosses correspond to the modes 
localized in the vicinity of the 6 -bounce PO and follow a 
well-distinguished behavior as well. The crosses in shades 
of blue correspond, as pointed out by the inset images la¬ 
belled from A to E, to quasi-periodic marginally stable 
orbits localized on the boundaries of the cavity. The light 
blue crosses correspond to modes with a single crown (A), 
as the medium and dark blue crosses correspond respec¬ 
tively to double ( B ) and triple ( D ) crowns. When the 
crosses intersect, the corresponding modes are a mix of 
single and double crown modes (resp. double and triple) 
as shown for the mode C (resp. E). It is important to 
note that this behavior of strong deviation from the stan¬ 
dard value IPR=3 is commonly encountered is systems 
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FIG. 7. Inverse Participation Ratio (IPR) as a function of the 
transverse wave number ft x a. (a) Chaotic case (7 = 3 / 2 ) ; 
As expected, the IPR concentrates in the vicinity of IPR=3 
except for the modes built on neutral orbits (upper inset, 
A). Inset B shows a generic ergodic mode, (b) Mixed case 
(7 = 0.95) ; Families of modes are distinguishable as shown by 
the different markers: red circles correspond to 2-PO regular 
modes, from red to green dots correspond to higher order reg¬ 
ular modes denoted from m = 1 to m — 9, purple crosses cor¬ 
respond to modes built on the 6-PO, light, medium and dark 
blue crosses correspond to modes localized on quasi-periodic 
marginally stable orbits (insets A to E). 


presenting a strong enough disorder to promote the ex¬ 
istence of (Anderson-) localized modes. Here the order , 
by means of the regular modes, is surprisingly responsi¬ 
ble for the deviation of the IPR. The modes correspond¬ 
ing to the lower IPRs - marked in black, right above 
IPR=3 (red dashed line) - correspond to chaotic modes 
that present an ergodic behavior. The density of points 
around IPR=3 is much less important than above, which 
confirms that ergodic modes are a minority in a mixed 
system characterized by the parameter 7 = 0.95. 

By use of a heuristic model, we derive an analytic ex¬ 
pression for the value of IPR for the high-ft regular modes 
along 2-PO with respect to fto, p . We assume that the 
modes are spatially localized on a surface <Si oc . Thus, the 
intensity of each mode equals Iq on <Si oc , and 0 elsewhere. 



FIG. 8. IPR as a function of ft x a for the regular modes 
and for high values of ft. The dots and circles correspond to 
the points presented in Fig. 0Jb) and the continuous lines 
correspond to the expression CG2>- 


Then, Iq = S/S\ oc and using equation (TUIl) . one gets: 

IPR = A- (12) 

<^loc 

The analogy of these modes with the fundamental gaus- 
sian modes uq :P of the stable resonator (see section ITTT|) . 
allows one to derive an analytic expression for S\ oc \ 



— w 0 L loc (13) 

where L\ oc is an effective “localization length”. Finally, 
one gets: 


IPR = 



(14) 


We obtain a perfect agreement between this expression 
and the IPR of the 2-PO regular modes ?i m =o,p [large red 
circles and red continuous line in Fig. [8]. 

Moreover, for large values of ft, the IPR appears to struc¬ 
ture itself by presenting a regular evolution for each fam¬ 
ily of modes. The IPR for the modes of a given m value 
are following a monotonous curve which is perfectly ad¬ 
justed by an empirical expression: 

IPR( / ft m? p) = c m n u (13) 

with = —0.017 xm+1 and c m a constant. 

In Figure [8j the red circles correspond to the 2-PO 
regular modes uo^ p and the line corresponds to the plot of 
the expression m ■ The dots in shades of colors from red 
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FIG. 9. Experimental setup. P: polarizer; M: mirror; A/2: 
half-wavelength waveplate; xlO, x20, x40: microscope ob¬ 
jectives; SLM: Spatial Light Modulator; L: lens needed to 
collect the FF; CCD: CCD camera. 


to green are the IPRs of higher order modes as shown in 
Fig. m and the lines are given by plotting the expression 
m- As one can note, the agreement between the values 
of the IPR and the analytic curves is very good. 


V. EXPERIMENTAL OBSERVATION AND 
MANIPULATION OF THE REGULAR MODES 
OF THE MIXED FIBER 

In order to study experimentally the light propagation 
into a mixed fiber, we manufactured a multimode op¬ 
tical fiber whose transverse shape is a truncated circle 
with 7 = 0.95. The illumination is made with a 2mW 
HeNe cw laser (Ao = 633nm). The polarization is con¬ 
trolled through a polarizer and a A /2 wave-plate. The 
laser beam, of initial diameter 2 mm, is extended in or¬ 
der to use all the surface of a Spatial Light Modulator 
(SLM, Amplitude only). The SLM is used to shape the 
intensity profile of the illumination beam. The beam is 
focused into the fiber thanks to a microscope objective 
in order to fit the surface of the fiber’s core. The angle 
of the initial beam is controlled by the angle of the input 
end of the fiber with respect to the optical axis. It allows 
to fix the average transverse wavenumber of the illumi¬ 
nation and thus to select the range of excited modes. A 
spatial modulation can be imposed through the SLM in 
order to shape the spatial distribution of the field inten¬ 
sity and then select a family of excited modes. At the 
output, the beam is collected by a microscope objective 
and sent to a CCD Camera. The NF is then directly col¬ 
lected. A supplementary lens is needed in order to collect 
the FF. 

At first, we did not use the SLM, and we only con¬ 
trol the illumination by the position of a focused beam 
in front of the input end of the fiber. Figures [TO^ a) and 
(g) show two typical excitations. The gaussian shape is 
obtained at the focal point of a microscope objective, and 
a tilt is given in order to select the range of transverse 
wave vectors propagating into the fiber. Figure [TOTbi 


presents the NF measured at the fiber’s output for the 
first illumination [Fig. ITUT a)]. The spatial repartition 
of the NF is characteristic of a superposition of regu¬ 
lar modes of the 2 -PO. One can see the large amount 
of intensity near the straight boundary which is a sig¬ 
nature of the excitation of these modes. Indeed, all the 
regular modes along the 2-PO present such a high inten¬ 
sity in the vicinity of the straight boundary, along the 
x-axis. Moreover, one can note the circular arc-like pat¬ 
tern also present on Figures [5][ a) and (e). Figure [TOT ci 
presents the FF intensity pattern of the outgoing wave. 
The anisotropic configuration along the k x -axis confirms 
that the field builds up on the 2-PO. The spatial ex¬ 
tension of the FF testifies that numerous gaussian-like 
modes of the type uo :P (x,y) have been excited. Figures 
ITOlf d) and (e) show the NF and FF resulting from nu¬ 
merical simulations. Using a standard beam propagation 
method algorithm [39j, we simulate the propagation of 
an initial gaussian beam along the optical fiber. The 
parameters here are given in order to reproduce the ex¬ 
perimental initial condition. As one can see, the agree¬ 
ment between the numerical and experimental results is 
very good, without any adjustable parameter. Figure 
ITOlf f) presents the numerical spectrum calculated through 
the standard method [40]. More precisely, \C(k)\ 2 is the 
square modulus of the Fourier transform of the correla¬ 
tion function defined as the overlap between the outgoing 
field and the initial condition. The observed Fabry-Perot- 
like spectrum confirms that the analogy between the 2 - 
PO regular modes of the mixed fiber and the gaussian 
modes of a stable optical cavity is relevant and one can 
measure a free spectral range S = (AA) xa = 7ra/d = 7r/7 
where Ak = K m ,p +1 — Km,p as defined in equation ([9|). 

Figure ITOT h) presents the NF collected for the second 
illumination [Fig. fTOf gi]. Here, the pattern is completely 
different and the spatial distribution of intensity is sta¬ 
tistically uniform, made of grains of light of random am¬ 
plitude and size. This speckle-like pattern results from 
the superposition of ergodic modes [see Fig. Eti)]. The 
ff [Fig. nurn] confirms this assumption, by presenting 
an isotropic distribution of wavevectors whose modulus 
are radially confined around an average wavenumber R. 
Figures flQlf i) and (k) present the associated numerical 
simulations and show a good agreement with the exper¬ 
imental results. The numerical spectrum [Figure fTUTl)] 
does not present any predominant structuration. More¬ 
over, it shows no discrimination of the modes: all the 
modes in the range of k corresponding to the illumina¬ 
tion are excited. The three regularly-spaced peaks one 
can see around k x a = 90 are a reminiscence of a periodic 
orbit that have been excited and which is responsible for 
the caustic visible on Figure fTUT i). 

The evaluation of the IPR from experiments for a 
superposition of modes is complicated due to the beat¬ 
ing between the excited modes along the propagation. 
However, it is always possible to get an estimation of 
the IPR at the fiber’s output and to compare the result 
with the numerical simulations for given propagation 
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FIG. 10. Initial illumination beams, NF and FF observed at the output of the fiber, corresponding numerical simulations and 
numerical calculated spectra, for an illumination favorable to the excitation of a family of (a-f) regular modes and (g-1) ergodic 
modes. 


lengths. We note that the IPR is always around 3 for 
the ergodic-like spatial repartition of the intensity and 
widely above for the superposition of regular modes of 
the 2-PO. This result confirms the very different nature 
of these two kinds of behaviors. 

From the spectra, we aim to extract some information 
on the underlying geometrical signature of the modes. 
The spectra presented in Figures [TOf fl and (1) are, as a 
matter of fact, equivalent to the density of states n(ft) so 
as n(^) = |C(ft)| 2 . An other expression of this density of 
states reads [4l| : 

n(ft) = no (ft) + ^^tjAj exp (inij) (16) 

j 

where the sum is running over the periodic orbits, ij is 
the total transverse optical length of orbit j and Aj en¬ 
compasses a classical amplitude related to the stability 
of the orbit and a phase associated to caustics and reflec¬ 
tions. The first part of the density, no (ft) is given by the 
so-called Weyl formula [42j : 


n 0 (ft) 


S_ _V_ 

27r 47 t 


(17) 


where V is the perimeter of the cross section of the fiber. 
As one can see in equation CG3), a relevant way to charac¬ 
terize the presence or not of any regularity in the under¬ 
lying ray dynamics is to calculate the “length spectra”, 
as the Fourier transform of n(ft) in order to get the quan¬ 
tity C(i) as a function of A C(t) will display peaks at the 


corresponding orbit length. Thus, the geometrical length 
of the periodic orbits can be directly extracted from the 
measure of the spectrum n(ft) which is calculated inde¬ 
pendently from the knowledge of the underlying periodic 
orbits (43| . 

Figure Era) presents the length spectrum calculated 
from Fig. [TOf fl. The first peak corresponds to the length 
of the 2-PO (pointed out by the red dashed line), that is 
f/a = 2 X7 = 1.9 (the factor 2 appearing for the round 
trip) and the other peaks being the harmonics. It is thus 
obvious that the only contribution to the spatial distribu¬ 
tion of the field is due to a superposition of regular modes 
of the 2-PO. Figure ITTfb) shows the length spectrum as¬ 
sociated to the spectrum of Fig. nnn). Some peaks still 
appear as residual resonances not predominant in the dy¬ 
namics. The length spectrum with the same scale as the 
latter is shown in the inset and it is clear that no specific 
trajectories are promoted. 

As shown right above, a beam focused close to the 
truncation - which corresponds to the location of the 
gaussian transverse modes waist - leads to the excitation 
of 2-PO regular modes. The width of the spectrum (the 
number of excited modes) then depends on the size of the 
illumination beam. A fine tuning of the illumination is 
obtained by adding a spatial modulation along the x-axis 
by use of a SLM: it allows one to select the order of the 
excited modes. Figures [T2T a and b) show examples of the 
FF for two different modulations associated to two differ¬ 
ent transverse wavenumbers ft. We notice that the mean 
radius of the intensity repartition in the far-field, that is 
the mean value of ft, is slightly different. A low-resolution 
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FIG. 11. Length spectra corresponding to spectra presented 
in (a) Fig. HEf): one can see the regularly distributed peaks 
corresponding for the first to the length of the 2-PO and for 
the others to the harmonics ; and in (b) Fig. I10IT ): one can see 
that no specific length exists (see inset) even if a zoom shows 
some peaks distributed in an arbitrary way. The arrows point 
the 2-PO and 5-PO respectively. 



FIG. 12. (a,b) Experimental FF measured for two different 
spatial modulations of the illumination beam, (c) Angular 
integrated FF. Blue line corresponds to case (a) and red line 
to case (b). 

tion, we slightly change the location of the maximum of 
intensity of the illumination beam and then excite pref¬ 
erentially regular modes with inverse parity. 

We perform a numerical simulation to analyze the ex¬ 
perimental angular far-field integration. We project two 
initial conditions analogous to the ones used in the exper¬ 
iment - that is an asymetric gaussian beam stretched and 
modulated in the ^-direction - on the calculated modes 
to get the weights c n of each mode in the propagating 
field y , z) describes as : 


spatial frequency spectrum is obtained from the angular 
integration of the far-field [Fig. H2lc)]. It allows to get 
rough informations on the transverse wavenumbers char¬ 
acterizing the propagated field. The blue line (resp. red 
line) corresponds to the FF presented in Fig. fl2T a) (resp. 

(b))- 

In Fig. [12j modulations clearly appear both in the /“in¬ 
direction in the FFs and in their angular integrations. It 
points out that only some specific modes (the 2-PO reg¬ 
ular modes) have been selected. By superimposing the 
expected values of k for the 2-PO regular modes calcu¬ 
lated for our experimental system using eq. m , one can 
notice that most of the maxima in the angular integrated 
far-field correspond to some k o )P values associated to the 
uo :P (x,y) 2-PO regular modes. As confirmed by the 2D 
representations of the FF [Figs. [T2l a,b)]. 2-PO regular 
modes are essentially excited. By changing the modula¬ 


fp(x, y,z) = ^2 Cn<t>n{x, y) exp (■ -jfi n z ) (18) 

n 

where 0 n is a guided mode and f3 n = y/k^ — its corre¬ 
sponding constant of propagation. Numerically, we can 
superimpose the angular integrated far-field with the cal¬ 
culated weights c n of each modes. Figure fl3T a) shows 
a behavior analogous to the experimental angular inte¬ 
grated far-fields. The peaks observed for low value of 
KXfl are also present and pointed by the vertical arrows. 
Actually, they also appear in the angular integrated far- 
field calculated for a unique 2-PO regular mode and are 
a consequence of the gaussian beam nature of the 2-PO 
regular modes along the stable trajectory. To compare, 
scar modes, which are the result of constructive interfer¬ 
ences between two counter-propagative plane waves, are 
characterized by two symmetric peaks in the FF. Con- 
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Hi x a 


FIG. 13. Angular integrated FF obtained by numerical simu- 
mations. (a) For an illumination favorable to the excitation 
of regular modes. Inset shows the FF ; (b) For a modula¬ 
tion along the ^-direction, so that no regular mode is being 
excited. Inset shows the FF. 


sequently, a unique peak would appear in the angular 
integrated FF. In Fig. [13j we observe a broaden peak in 
the integrated FF, which is due to the interferences be¬ 
tween the excited modes preferentially being 2-PO regu¬ 
lar modes. Finally, we report in Fig. [T3l bk the angular 
integrated far-fleld for an initial modulated beam out of 
the 2-PO direction (precisely, in the upper part of the 
fiber, with a modulation having the same k value but 
along the ^/-direction to avoid the 2-PO). We observe the 
vanishment of the first peaks associated to the character¬ 
istic signature of the 2-PO regular modes and an enlarged 
peak that corresponds to a superposition of a great num¬ 
ber of arbitrary excited modes around the main n value. 


VI. CONCLUSION AND DISCUSSION 

In this paper we have presented numerical and experi¬ 
mental investigations of the modes of a highly multimode 
fiber whose transverse cross section is designed to induced 
a complex ray dynamics. This multimode optical fiber 
is clearly an ideal experimental support for this kind of 
study, as it allows an easy-to-implement setup, with a 
controllable excitation and a direct visualization. More¬ 
over, this system is versatile in the sense that a slight 
change in the level of the truncation leads to the explo¬ 
ration of different types of dynamics. One thus has ac¬ 
cess, with the same experimental system, to mixed, reg¬ 
ular as well as chaotic dynamics. Here, we focused our 
work on the study of a mixed dynamics. In particular, we 
concentrated the analysis on a family of regular modes of 
the system, that presents a spatially localized intensity 
pattern, as well as a Fabry-Perot-like discrete spectrum. 
An analogy with the gaussian modes of an optical res¬ 
onator leads to the derivation of an analytical expression 
for the regular modes as well as a resonant condition 
justifying the regular distributed spectrum. The strong 
spatial localization of the regular modes finds some sig¬ 
natures in the study of the Inverse Participation Ratio, 
by largely exceeding the value predicted by the Random 
Matrix Theory. By coupling these results with the anal¬ 
ogy with the gaussian modes, we then were able to give 
an expression for the IPR as a function of the wavenum¬ 
ber for the regular modes. 

With a suitable shaping of the initial beam, we 
experimentally demonstrated that these modes can be 
selectively excited, and that they are robust to mode 
coupling. This ease of shielding the other modes open 
the way to multiple applications in optical telecommu¬ 
nications. For instance, the regular modes of the mixed 
optical fiber are perfectly suitable for Mode Division 
Multiplexing (MDM) and would allow the use of much 
more transmission channels than already achieved in 
conventional devices. Moreover, these modes can benefit 
from a selective optical amplification, by optimizing 
the spatial overlap with a gain medium. To do so, one 
simply has to locate the active medium in the vicinity 
of the truncation, where the regular modes have their 
maximum of intensity 0 . In a more fundamental 
point of view, a nonlinear mixed optical fiber would 
allow an enhancement of the phenomenon of optical 
thermalization and condensation of classical waves, as 
the spatial overlap between the modes is an important 
parameter ji3 |. 

The authors thank gratefully M. Ude and S. Trzesien 
for their involvement in the preform and fiber manufac¬ 
ture. 
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